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Abstract 

The adsorption of highly oppositely charged flexible polyelectrolytes (PEs) on a charged planar 
substrate is investigated by means of Monte Carlo (MC) simulations. We study in detail the equi- 
librium structure of the first few PE layers. The influence of the chain length and of a (extra) 
non-electrostatic short range attraction between the polycations and the negatively charged sub- 
strate is considered. We show that the stability as well as the microstructure of the PE layers are 
especially sensitive to the strength of this latter interaction. Qualitative agreement is reached with 
some recent experiments. 
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I. INTRODUCTION 



PE multilayer structures are often obtained in a so-called layer-by-layer method by al- 
ternating exposure of a charged substrate to solutions of polycations (PCs) and polyanions 
(PAs). This widely used technique was first introduced by Decher and coworkers and 
its simplicity and versatility trigger a large interest in the engineering community As exam- 
ples of technological applications, one can mention: biosensing [3], catalysis [4], non-linear 
optical devices 3] , nanoparticles coating 0, Q| , etc. 

On the theoretical side, there exit a few analytical works about PE multilayers on charged 
planar surfaces based on different levels of approximation (8,0, n|. Solis and Olvera de la 
Cruz considered the conditions under which the spontaneous formation of polyelectrolyte 
layered structures can be induced by a charged wall [8]. Based on Debye-Hiickel approx- 
imations for the electrostatic interactions and including some lateral correlations by the 
consideration of given adsorbed PE structures, Netz and Joanny Q] found a remarkable sta- 
bility of the (semi-flexible) PE multilayers supported by scaling laws. For weakly charged 
flexible polyelectrolytes at high ionic strength qualitative agreements between theory |10(, 
also based on scaling laws, and experimental observations [llj (such as the predicted thick- 
ness and net charge of the PE multilayer) have been provided. 

The important driving force for all these PE multilayering processes is of electrostatic 
origin. More precisely, it is based on an overcharging mechanism, where the first layer over- 
charges the substrate and, along the PE multilayering process, the top layer overcharges 
the adsorbed PE layers underneath. Nevertheless, the strong correlations existing between 
oppositely charged polyions, especially for highly charged PEs, provide a formidable chal- 
lenge for the understanding of the PE multilayer microstructures. In this respect, numerical 
simulations are of great help. It is only recently, that MC simulations were carried out to 
study such PE structures built up on spherical charged surfaces [12j . 

In this paper, we provide a detailed study of the PE multilayer structure adsorbed on a 
charged planar surface and discuss the basic mechanisms that are involved there by means 
of MC simulations. Our paper is organized as follows: Sec. HH is devoted to the description 
of our MC simulation technique. The measured quantities are specified in Sec. IIHI The PE 
monolayering is studied in Sec. IIVI and the PE bilayering in Sec. |V] Then the PE multi- 
layering process is addressed in Sec. IVI1 Finally, Sec. I VI II contains some brief concluding 
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remarks. 



II. SIMULATION DETAILS 

The setup of the system under consideration is similar to that recently investigated with 
a spherical substrate Within the framework of the primitive model we consider a 

PE solution near a charged hard wall with an implicit solvent (water) of relative dielectric 
permittivity e r ~ 80. This charged substrate located at z = is characterized by a negative 
surface bare charge density — cr e, where e is the (positive) elementary charge and a > 
is the number of charges per unit area. Electroneutrality is always ensured by the presence 
of explicit monovalent (Z c = 1) plate's counterions of diameter a. PE chains (iV + PCs and 
AL PAs) are made up of N m monovalent monomers (Z m = 1) of diameter a. Hence, all 
micro ions are monovalent: Z = Z c = Z m = 1. For the sake of simplicity, we only consider 
here symmetrical complexes where PC and PA chains have the same length and carry the 
same charge in absolute value. 

All these particles making up the system are confined in a L x L x r box. Periodic 
conditions are applied in the (x, y) directions, whereas hard walls are present at z = 
(location of the charged plate) and z = r (location of an uncharged wall). To avoid the 



appearance of image charges jjl3 



14j , we assume that on both parts of the charged plate (at 



z — 0) the dielectric constants are the same. 

The total energy of interaction of the system can be written as 



Utot = £ K ate) (^) + U^ e \zi) + U^T\ Zl )\ + (1) 

i 

E Phsinj) + U cou i(r i:j ) + U FE NE(rij) + U L j(r i:j )} , 

i,i<j 

where the first (single) sum stems from the interaction between an ion % (located at z = zi) 
and the charged plate, and the second (double) sum stems from the pair interaction between 
ions i and j with = |r, ; — rj\. All these contributions to U tot in Eq. (JTJ) are described in 
detail below. 
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Excluded volume interactions are modeled via a hardcore potential defined as follows 



oo, for Tij < a 



0, for Tij > a 



(2) 



for the microion-microion one, except for the monomer- monomer one |15|], and 



U, 



(plate) 



1 1 s 



Zi) 



oo, for Zi < a/2 

oo, for Zi > r — a/2 

0, for a/2<Zi<r- a/2 



(3) 



for the plate-microion one. 

The electrostatic energy of interaction between two ions i and j reads 



U C oul yij / 



±- 



U 



(4) 

where +(-) applies to charges of the same (opposite) sign, and Ib = e 2 / '47reoe r /csT is the 
Bjerrum length corresponding to the distance at which two monovalent ions interact with 
ksT. The electrostatic energy of interaction between an ion i and the (uniformly) charged 
plate reads 

c t J: i} = ±2nl B a oZi (5) 

Kb-L 

where +(-) applies to positively (negatively) charged ions. An appropriate and efficient 
modified Lekner sum was utilized to compute the electrostatic interactions with periodicity 



in two directions 



16l Il7| . To link our simulation parameters to experimental units and 



room temperature (T = 298K) we choose a = 4.25 A leading to the Bjerrum length of water 
Ib = 1.68a = 7.14 A. The surface charge density of the planar macroion was chosen as 
-a e w -0.165 C/m 2 . 

The polyelectrolyte chain connectivity is modeled by employing a standard FENE poten- 
tial in good solvent (see, e.g., [18]), which reads 



UfeneIt) 



-~kRI Id 



Rl 



oo, 



for r < Rq 



for r > Rn 



(6) 



with k = 27ksT/a 2 and Rq = 1.5a. The excluded volume interaction between chain 
monomers is taken into account via a purely repulsive Lennard- Jones (LJ) potential given 



by 



U LJ {r) 



4e 



T. 



+ e, for r < 2 1 / 6 c 



(7) 



0, for r > 2 1 / 6 a 

where e = /cgT. These parameter values lead to an equilibrium bond length I = 0.98a. 

An important interaction in PE multilayering is the non-electrostatic short ranged at- 
traction, U^ te \ between the planar macroion and the PC chain. To include this kind 
of interaction, we choose without loss of generality a (microscopic) van der Waals (VDW) 
potential of interaction between the planar macroion and a PC monomer that is given by 



U. 



(plate) 



vdw 



for z > a/2, 



(8) 



z + a/2 

where Xvdw is a positive dimensionless parameter describing the strength of this attraction. 
Thereby, at contact (i.e., z = a/2), the magnitude of the attraction is Xvdwt = XvdwksT 
which is, in fact, the relevant characteristic of this potential. Since it is not straightforward 
to directly link this strength of adsorption to experimental values, we choose Xvdw = 5 (also 



12j), so as to mimic good 



considered among other values in the case of a spherical macroion 
"anchoring" properties to the planar substrate. 

All the simulation parameters are gathered in Table The set of simulated systems can 
be found in Table ITT1 The equilibrium properties of our model system were obtained by using 
standard canonical MC simulations following the Metropolis scheme 0,13- Single-particle 
moves were considered with an acceptance ratio of 30% for the monomers and 50% for the 
counterions. Typically, about 5 x 10 4 to x 10 6 MC steps per particle were required for equili- 
bration, and about 5 x 10 5 — 10 6 subsequent MC steps were used to perform measurements. 
To improve the computational efficiency, we omitted the presence of PE counterions when 
iV + = i\L so that the system is still globally electroneutral. We have systematically checked 
for iV + = iV_ = 20 (system C) that the (average) PE configurations (especially the monomer 
distribution) are indistinguishable, within the statistical uncertainty, from those where PE 
counterions are explicitly taken into account, as it should be. 
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III. MEASURED QUANTITIES 

We briefly describe the different observables that are going to be measured. In order to 
characterize the PE adsorption, we compute the monomer density n±(z) that is normalized 
as follows 

pr-a/2 

I n±(z)L 2 dz = N±N m (9) 

J a/2 

where +(— ) applies to PCs (PAs). This quantity is of special interest to characterize the 
degree of ordering in the vicinity of the planar macroion surface. 

The total number of accumulated monomers N±(r) within a distance z from the planar 
macroion is given by 

N±(z) = [ n±(z')L 2 dz' (10) 

J a/2 

where +(— ) applies to PCs (PAs). This observable will be addressed in the study of PE 
monolayer (Sec. HVjl and PE bilayer (Sec. Ivj) . 

Another relevant quantity is the global net fluid charge a(z) which is defined as follows 

a(z)= [ [h+{z')-n„{z')]dz', (11) 

J a/2 

where h + (fi_) stand for the density of all the positive (negative) microions (i.e., monomers 
and counterions) . Thus, a(z) corresponds to the net fluid charge per unit area (omitting 
the bare macroion surface- charge — er ) within a distance z from the charged wall. At the 
uncharged wall, electroneutrality imposes a(z = r — a/2) = a~o- By simple application of 
the Gauss' law, [u{z) — <jq\ is directly proportional to the mean electric field at z. Therefore 
a(z) can measure the screening strength of the macroion-plate charge by the neighboring 
solute charged species. 

IV. MONOLAYER 

In this part, we study the adsorption of PC chains for two chain lengths N m = 10 (system 
A) and N m = 20 (system B), and for two different couplings Xvdw = and Xvdw = 5. 
Experimentally, this would correspond to the formation of the first polyelectrolyte layer. 
This is a decisive step to elucidate the even more complex PE multilayer structures where 
additionally PAs are also present. 
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Here, where AL = (i.e., no polyanions), global electroneutrality is ensured by the 
presence of explicit PC's counterions (i.e., monovalent anions) and the macro ion-plate's 
counterions (i.e., monovalent cations). Also, we recall that the total number of monomers, 
N + N m = 200, is identical for both systems A and B under consideration (see Table In]) . 
Hence, the total monomer charge is the same for systems A and B. 

The profiles of the monomer density n+(z) are depicted in Fig. UJ At Xvdw = 0, the density 
n+(z) near contact (z ~ a/2) is basically independent on the chain size N m . But away from 
the surface, the density of monomers is slightly higher for larger N m . This is a combined 
effect of (i) entropy and (ii) electrostatic correlations. These underlying mechanisms at 
Xvdw = can be explained with simple ideas as follows: 

• At fixed number of total monomers, entropic effects are larger the shorter the chains, 
and in the limiting case of N m = 1 (i.e., the electrolyte limit) entropy effects are maxi- 
mal leading to the highest monomer "release" . It is to say that the chain connectivity 
lowers the entropy of the system. 



In parallel, electrostatic correlations 2l|, 2^, l^j are also higher the higher the valence 



of the adsorbed particles. In our case N m plays the role of the polyion valence. 

The density n+(z) near contact increases considerably with Xvdw (here about one order of 
magnitude) as expected. It turns out that with Xvdw — 5, the n + (^-profiles are basically 
identical for N m = 10 and N m = 20. This is due to the sufficiently strong non-electrostatic 
attractive force that can overcompensate the antagonistic entropic effects that were more 
efficient at Xvdw = 0. 

The fraction N + (z)/(N + N m ) of adsorbed monomers can be found in Fig. [2J At a 
z-distance of 1.5a from the planar macroion surface (corresponding to a width of two 
monomers), about 90% of the monomers are adsorbed for Xvdw — 5 against only ~ 50% 
for Xvdw = 0. Again, at given Xvdw, N+(z) / (N + N m ) is larger for longer chains due to the 
same coupled effects of entropy and electrostatic correlations explained above. 

The (global) net fluid charge o~(z) is reported in Fig. |H1 In all cases we observe a 
macroion-surface charge reversal (i.e., cr(z)/a > 1). The position z = z* at which <j{z*) 
gets its maximal value decreases with Xvdw, due to the Xudu>- ermance d adsorption of the 
PCs. Concomitantly, this overcharging increases with Xvdw, since the (extra) gain in energy 
by macro ion- monomer VDW interactions can better overcome (the higher Xvdw) the cost 



of the self-energy stemming from the adsorbed excess charge |l2j. More quantitatively, we 
have a(z*)/a « 1.7 at Xvdw = 5 against only a(z*)/o~o ~ 1.25 at Xvdw = 0. Note that the 
maximal value of charge reversal of (200 — 90)/90 = 122% (i.e., a(z*)/a = 2.22) allowed 
by the total charge of PCs can not be reached due to a slight accumulation of microanions. 
In agreement with the profiles of n+(z) and N+(z) (see Fig. [T]and Fig. |2J), at given Xvdw, 
the overcharging gets higher the higher the chain length. Those (locally) overcharged states 
should be the driving force for the building of subsequent PE bilayers when PA chains are 
added. 

Typical equilibrium configurations can be found in Fig. The qualitative difference 
between Xvdw — [Figs. |Ua) and (b)] and Xvdw = 5 [Figs. EJc) and (d)] is rather spectacular. 
Without additional VDW attraction (x V dw = 0) the adsorption is much weaker than at 
Xvdw = 5, where in the latter situation the ^-fluctuation is very weak within the adsorbed 
layer. Basically the first layer is glued at Xvdw — 5, and the excess PC chains float in 
the solution. It is typically this type of configurations for the first layer that is wanted in 
experimental situations. 

The next section (Sec. E} that concerns bilayering will show that the (enhanced) stability 
of this first layer is decisive for the onset of multilayers. 

V. BILAYER 

We now consider the case where additionally PA chains are present (systems C and D), 
so that we have a neutral polyelectrolyte complex (i.e., N + N m = N_N m = 200). Global 
electroneutrality is ensured by the counterions of the planar macroion as usual. For such 
parameters, the final equilibrium structure consists essentially of bilayers with sometimes 
the onset of a weakly stable third layer. Experimentally this would correspond to the process 
of the second polyelectrolyte layer formation (with system A or B as the initial state). We 
stress the fact that this process is fully reversible for the parameters investigated in our 
present study. In particular, we checked that the same final equilibrium configuration is 
obtained either by (i) starting from system A or B and then adding PAs or (ii) starting 
directly with the mixture of oppositely charged polyelectrolytes. 

The profiles of the monomer density n±(z) at Xvdw = and Xvdw = 5 are depicted in 
Fig. Efa) and Fig. Efb), respectively. The corresponding microstructures are sketched in 
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Fig. El The density of PC monomers n + (r) near contact increases considerably with Xvdw 
as expected. Interestingly, at Xvdw = 0, a comparison with systems A and B (see Fig. |TJ) 
indicates that the adsorption of PC monomers is weaker when additional PAs are present. 



This effect was already observed with spherical substrates 



12 1, and the same mechanisms 



apply here to planar surfaces. More explicitly, the PC chain tends to build up a globular state 
(reminiscent of the classical bulk PE collapse j^J) by getting complexed to the PA chain, as 
well illustrated in Fig. HJa) and Fig. Efb). Thereby, the mean monomer coordination number 
(or the mean number of monomer neighbors) gets higher which is both (i) entropically and (ii) 
energetically (at least from the PE complex viewpoint) favorable. This PC desorption is only 
appreciable at sufficiently low Xvdw where the energy loss stemming from the PC desorption 
is well balanced (or even overcompensated depending generally on the parameters) by the 
energy gained in building a PC-PA globular structure. This "auto-globalization" is also 
enhanced by increasing N m as it should be [compare Fig. Efa) and (b)]. Note also that 
there is a small second peak in n+(z) at z pa 3.8a [see Fig. Efa)], which is rather the 
signature of a strong PC-PA globalization than a third PE layer. Besides, the peak in the 
PA density n-(z) located at z = z* pa 2.3a [see Fig. ED^a)], which is relatively far from that of 
a compact bilayer where z* = 1.5a, indicates the diffuse character of the bilayer at Xvdw — 0. 

At Xvdw = 5, the n+ (^-profiles are basically identical for N m = 10 and N m = 20. In 
contrast to Xvdw — 0, n + (z) near contact is somewhat larger at Xvdw = 5 and it is going to 
be explained later in the discussion of N + (z). As far as the PA density ri-(z) is concerned, 
we see that the peak is roughly 2-3 times higher (depending on N m ) with Xvdw = 5 than 
with Xvdw = 0. Also, its position (z* pa 1.5a) corresponds to that of a compact bilayer. A 
visual inspection of Fig. (HJ^c) and (d) confirms this feature. This again shows how important 
is the role of extra non-electrostatic attractive force for the stability of bilayers. 

An intermediate conclusion can be drawn from the above findings and especially from 
the microstructures depicted in Fig. |HJ 

• True bilayering (i.e., flat and dense layers) can only occur at non-zero Xvdw, as already 



reported for spherical charged substrates 



12j with large curvature. 



An interesting common characteristic of the microstructures at Xvdw — and Xvdw = 5 is 
the formation of small islands (along the substrate) made up of more or less flat (depending 
on Xvdw) PC- PA complexes, easily identifiable at N m = 20 [see Fig. MJd) and Fig. EJd)]. 
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The fraction N+(z) / \N + N m ) of adsorbed monomers at Xvdw = and Xvdw = 5 can be 
found in Fig. [7(a) and (b), respectively. A close look at Fig. UJa) reveals a smaller PC 
monomer accumulation (at Xvdw = 0) up to z ~ 3a (independently of N m ) than in the 
case where PA chains were absent (compare with Fig. |2J). This is fully consistent with the 
formation of PC-PA globules (relevant at Xvdw = 0) leading to the effective PC desorption 
already discussed above. In parallel, this PC-PA globalization tends to cancel the effect of 
chain length N m on N + (z). On the other hand, at Xvdw = 5, the situation is qualitatively 
different where the presence of PAs now induces an increase of N+(z) [compare Fig. UJb) 
and Fig. |2j. This phenomenon can be explained by electrostatic correlation effects. Indeed, 
at Xvdw — 5, the highly stable PC layer attracts more PA monomers than at Xvdw — 0, and 
thereby, "super" dipoles made of PC- PA monomer pairs build up, that are perceptible in 
Fig. Uc) and Fig. E(d). This leads to a strong attractive correlation interaction between 
the plate and those dipoles. In other terms the effect of finite Xvdw is to (strongly) polarize 
the adsorbed charged chains. Note also that a Xvdw = 5 the PC-PA globalization is much 
less favorable than at Xvdw = due to the higher cost of PC desorption energy in the former 
case. As a net effect there can be more adsorbed PC monomers compared to Xvdw = 0. In 
that case of Xvdw — 5, it is precisely this mechanism that tends to cancel the effect of N m 
on N + (z). As far as the PA monomer fraction N-(z) is concerned, Fig. UJ shows that the 
adsorption of monomers is much weaker and more diffuse at Xvdw = than at Xvdw = 5, as 
expected from Figs. El and El 

The net fluid charge a(z) is reported in Fig. El In all cases, the planar macroion gets 
overcharged and undercharged as one gets away from its surface. That is we have to deal 
with charge oscillations. Our results clearly show that the amplitude of those oscillations 
is systematically larger at high N m , as also observed without PAs (see Fig. |3J). This is 
consistent with the idea that lateral electrostatic correlations are enhanced by increasing 
the valence of the polyions (here N m ). Nevertheless, as soon as oppositely charged polyions 
can interact, there is a subtle interplay between clustering and the lateral correlations of 
polyions that governs the degree of overcharging near the planar macroion. At Xvdw = 5, we 
observe a significantly higher overcharging than without PAs (compare with Fig. |3J). This is 
in agreement with the profiles of N + (z) discussed previously. However, the positions of the 
first peak (z* ~ a for Xvdw — 5 and z* « 1.8a for Xvdw — 0) in a(z) remain nearly unchanged 
by the presence of PAs (compare Fig. [HI with Fig. 
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VI. MULTILAYER 



Presently, we consider the case where there are enough polyelectrolytes (N + N m = 
N_N m = 400) in the system to produce multilayers (systems E and F). Hence, com- 
pared to systems C and D, we have now doubled the polyelectrolyte concentration. Global 
electroneutrality is ensured by the counterions of the planar macroion as usual. 

The density profiles of n±(r) for Xvdw — and Xvdw = 5 are depicted in Fig. |§fa) and Fig. 
EJb), respectively. The corresponding microstructures are sketched in Fig. [HI In general, 
the densities of PC and PA monomers are systematically larger than those found for systems 
C and D corresponding to a lower PE concentration. This effect is due to the fact that, 
at higher concentration of oppositely charged chains, the number of dipoles (i.e., PC-PA 
monomer pairs) are also larger, and from this it results larger plate-dipole correlations. 

Even at Xvdw = with N m = 10, we can observe a non- negligible second peak in n + (z) 
(located at z m 3.8a) which is the signature of a third layer. This finding contrasts with 
what was observed at spherical substrates |l2j (also with Xvdw = 0, N m = 10, and with a 
similar macroion surface charge density), where not even a stable bilayer could build up. 
This radically different behavior can be accounted by geometrical arguments. Indeed, the 
potential of electrostatic interaction scales like 1/r in spherical geom etry against z in planar 
one. Hence, at sufficiently high curvature (as it was the case in Ref. |l2j where N m a/r > 1 
(25! with ro being the radius of the spherical macroion), qualitative differences from the 
planar case are then expected. However, the corresponding microstructure [see Fig. HOT a)] 
suggests a relatively large formation of PC-PA globules leading to a diffuse and porous 
multilayer. Always at Xvdw = 0, but with longer chains (N m = 20), Fig. EJa) shows that 
the degree of layering is higher as expected. This feature is well illustrated by Fig. ITOT b). 
where the PA monomers are visibly more attracted to the planar macroion surface. 

At Xvdw = 5, the adsorption of monomers is drastically increased due to the enhanced 
stability of the first PC layer that, in turn, induces a larger adsorption of the subsequent PAs 
and PCs. Compared to Xvdw = 0, all the peaks in n±(z) are shifted to smaller z, indicating 
a higher compaction. These higher ordering and compaction at Xvdw = 5 can be visually 
checked in Fig. flOT c) and (d). 

The net fluid charge a(z) is reported in Fig. ^2 As expected charge oscillations are 
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detected. However in all cases, the corresponding amplitude is decaying (as also found in 



Ref. 12j for spherical geometry). Compared to the bilayer situation (see Fig. |SJ), one 
remarks that the charge oscillations are now larger due to the enhanced "plate-dipole" 
correlations occurring at higher chain concentrations (as discussed above). On the other 
hand the positions of the extrema in the charge oscillations remain quasi unchanged. 

VII. CONCLUDING REMARKS 

We first would like to briefly discuss our findings with some experimental examples. 
Our results concerning the first layer (i.e., single PC layer) show that an additional non- 
electrostatic force is needed to enhance its stability. Experimentally, this is achieved by 
choosing PCs with good "anchoring" properties to a given substrate. In our model this was 
done by taking Xvdw — 5. This being said, the case Xvdw = is from a fundamental point of 
view interesting, since it corresponds to a purely electrostatic regime. 

Recently, Menchaca et al. found, by means of "liquid-cell AFM", that PE-complex grains 



appear at the first PE-layers [26|. This kind of structure (that we referred to as small islands 
- see Fig. |HJ) are confirmed by our simulations. Concomitantly, a significant roughness of 
the deposited bilayer was also detected in this experiment, which is directly linked to the 
presence of those grains. This microstructure seems also to be (indirectly) reported in other 
experiments using ellipsometry j^, where it is found that the structure of the two first 
bilayers are more porous than that of later bilayers. This is also in qualitative agreement 
with our microstructures depicted in Fig. |B] and Fig. However, more simulation data 
are needed to understand the PE structure beyond two bilayers. 

The degree of charge inversion of the substrate can be indirectly obtained by measuring 



the ^-potential via electrophore sis 



llj . In their experiment, Ladam et al. [11] observed that, 
the ^-potential profile is symmetrically oscillating. This 



after a few deposited PE layers [5 
reveals a "stationary" regime where, successively, polycations and polyanions are adsorbed 
with the same strength. Unfortunately, it is not possible for us to investigate numerically this 
regime due to the highly prohibitive computation time required there. However, the charge 
oscillations observed in our a(2;)-profiles indicate that by increasing the amount of layers, 
one first increases the amplitude of these oscillations. This confirms at least the general 
experimental evidence of the non-stationary regime a the early stage of PE multilayering. 
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We also would like to mention the possible effect of image charges stemming from the 
dielectric discontinuity between the substrate (typically e r « 2 — 5) and the solvent (here 
e r ~ 80), as is the case under experimental conditions. It is expected that image forces 
become especially relevant for PE monolayering (i.e., when PCs solely are present) 
Indeed for multivalent ions, a strong self-image repulsion occurs and leads to a shifted 
density-profile n + (z) with a maximum located somewhat further than the contact region 
[lllQl- However for PE multilayering, due to the presence of oppositely charged PEs the 
effect of image forces is considerably reduced (especially sufficiently away from the wall) due 
to the (self-)screening of the image charges. 

In summary, we have investigated by means of extensive MC simulations the equilibrium 
buildup of the few first layers adsorbed on a charged planar substrate. Two parameters were 
considered: (i) the chain length N m and (ii) the extra non-electrostatic short-range attraction 
(characterized here by Xvdw) between the planar macro ion surface and the polycation chains. 

For the bilayering, it was demonstrated that, within the electrostatic regime (i.e., 
Xvdw = 0), significant PC- PA globules build up leading to a very "porous" and diffuse 
bilayer structure. The PC- PA globalization is enhanced with N m . At sufficiently large Xvdw 
(here Xvdw — 5), the bilayer is much less diffuse and the oppositely charged chains are more 
polarized, leading to a high stability of the structure. 

The same qualitatively applies to the case of the two-bilayer (i.e., four PE layers) ad- 
sorption. Within this regime of layering as investigated here (up to four layers), we also 
found a non-linear regime, where for instance the separation of the peaks in the monomer 
densities are not identical. This is in qualitative agreement with the finding of Ladam et 



al. 



Ill ] where they reported a non-linear regime in the so-called "region I" corresponding 



to the PE multilayer-region close to the buffer |28j . The effect of N m is to globally enhance 
the stability of the multilayer structure due to the higher electrostatic correlations and also 
due to entropic effects. 

A future study should take into account the rigidity of the chain, which can drastically 
change the multilayer structure depending the stiffness. The formation of PE multilayers 
on cylindrical substrates seems also to be a promising research area, and to our knowledge 
it has never been investigated so far j^j. 
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TABLES 



TABLE I: Model simulation parameters with some fixed values. 



Parameters 



T = 298K room temperature 

(To = 90/L 2 macroion surface charge 

Z = 1 microion valence 

a = 4.25 A microion diameter 

Ib = 1.68a = 7.14 A Bjerrum length 

L = 22a (x, y)-box length 

r = 75a z-box length 

N + number of PCs 

number of PAs 

Npe = N + + iV_ total number of PEs 

N m number of monomers per chain 

Xvdw strength of the specific VDW attraction 



TABLE II: System parameters. The number of counterions (cations and anions) ensuring the 
overall electroneutrality of the system is not indicated. 



System N PE N + N_ N m 

A 20 20 10 

B 10 10 20 

C 40 20 20 10 

D 20 10 10 20 

E 80 40 40 10 

F 40 20 20 20 
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FIGURE CAPTIONS 

1. Profiles of PC monomer-density n+(z) at different Xvdw couplings (systems A and B). 
The inset corresponds to Xvdw = 5 where the two curves (N m = 10 and 20) are nearly 
indistinguishable. 

2. Fraction N + (z) / (N + N m ) of adsorbed PC monomers at different Xvdw couplings (sys- 
tems A and B). 

3. Net fluid charge a(z) at different Xvdw couplings (systems A and B). 

4. Typical equilibrium configurations for PC chains adsorbed onto an oppositely charged 
planar macroion (systems A and B). The little counterions are omitted for clarity. 

(a) Xvdw = 0, N m = 10 (b) Xvdw = 0, N m = 20 (c) Xvdw = 5, N m = 10 (d) Xvdw = 5, 
N m = 20. 

5. Profiles of monomer density n±(z) for oppositely charged polyelectrolytes (systems C 
and D). (a) Xvdw = 0. (b) Xvdw = 5. 

6. Typical equilibrium configurations for the adsorption of oppositely charged PE chains 
(systems C and D) onto a planar macroion. The polycations are in white and the 
polyanions in red. The little ions are omitted for clarity, (a) Xvdw = 0, N m = 10 (b) 

Xvdw = 0, N m = 20 (c) Xvdw = 5, N m = 10 (d) Xvdw = 5, N m = 20. 

7. Fraction of adsorbed monomers N±(z) for oppositely charged polyelectrolytes (systems 
C and D). (a) Xvdw = 0. (b) Xvdw = 5. 

8. Net fluid charge a(z) at different Xvdw couplings (systems C and D). 

9. Profiles of monomer density n±(z) for oppositely charged polyelectrolytes (systems E 
and F). (a) Xvdw = 0. (b) Xvdw = 5. 

10. Typical equilibrium configurations for the adsorption of oppositely charged PE chains 
(systems E and F) onto a planar macroion. The polycations are in white and the 
polyanions in red. The little ions are omitted for clarity, (a) Xvdw = 0, N m = 10 (b) 

Xvdw = 0, N m = 20 (c) Xvdw = 5, N m = 10 (d) Xvdw = 5, N m = 20. 

11. Net fluid charge a(z) at different Xvdw couplings (systems E and F). 
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FIG. 3: Messina, Macromolecules 
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FIG. 7: Messina, Macromolecules 
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FIG. 9: Messina, Macromolecules 
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FIG. 11: Messina, Macromolecules 
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